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(57) ABSTRACT 

Method and system for determining an optimal route for an 
aircraft moving between first and second waypoints in a 
general wind environment A selected first wind environ- 
ment is analyzed for which a nominal solution can be 
determined. A second wind environment is then incorpo- 
rated; and a neighboring optimal control (NOC) analysis is 
performed to estimate an optimal route for the second wind 
environment. In particular examples with flight distances of 
2500 and 6000 nautical miles in the presence of constant or 
piecewise linearly varying winds, the difference in flight 
time between a nominal solution and an optimal solution is 
3.4 to 5 percent Constant or variable winds and aircraft 
speeds can be used. Updated second wind environment 
information can be provided and used to obtain an updated 
optimal route. 

32 Claims, 4 Drawing Sheets 
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NEIGHBORING OPTIMAL AIRCRAFT 
GUIDANCE IN A GENERAL WIND 
ENVIRONMENT 

ORIGIN OF THE INVENTION 

The invention disclosed herein was made by employees of 
the United States Government and may be manufactured and 
used by or for the Government for governmental purposes 
without payment of any royalties for such manufacture and 
use. 

FIELD OF THE INVENTION 

This invention relates to optimization of a flight route of 
an aircraft in an environment with constant or variable 
winds. 

BACKGROUND OF THE INVENTION 

Commercial and private aircraft travel great distances, 
500 to 5,000 miles, at varying altitudes and are subjected to 2 o 
winds up to 100 knots (120 miles per hour) in varying 
directions. Over the course of a long flight, requiring several 
hours to complete, the aircraft may change course and/or 
change altitude to avoid a heavy storm or a second aircraft 
that is climbing through the first aircraft’s present cruise 25 
altitude. Even if the aircraft cruises at a fixed altitude, the 
prevailing winds will vary strongly in magnitude and direc- 
tion from one along the route. For typical wind patterns 
encountered across the North American continent (e.g., in a 
flight from New York to San Francisco), traveling a wind- 
optimized flight route can reduce flight time, relative to 
flight time over a great circle route, by up to 8 percent. For 
a flight from the JFK Airport in New York to the San 
Francisco Airport, this reduction in time can be 20 minutes 
or more. For shorter flights, the aircraft flight speed, the 
wind speed and wind shear are also lower and the relative 
effect of variable local winds is stronger, and the ratio of 
time required to travel a great circle route and time required 
to travel an optimized route is approximately constant. 

However, computation from first principles of an optimal 
flight route between first and second waypoints in the 
presence of variable, even piecewise linear, wind direction 
and speed requires complex calculations involving many 
variables and may require too much time to be practical. 
What is needed is a less complex approach that provides a 
near-optimal flight route and that can compensate for the 
changing wind conditions encountered by an aircraft in 
flight. Preferably, the approach should allow use of raw 
numerical data or of analytical approximations for wind 
speed and direction encountered along a flight path and 
should allow for variable aircraft speed and altitude along 
the flight path. 

SUMMARY OF THE INVENTION 

These needs are met by the invention, which begins with 
a nominal solution that satisfies the source and destination 
waypoint constraints and applies a neighboring optimal 
control (“NOC”) and guidance analysis to a wind velocity 
vector field and aircraft speed prescription and provides a 
near-optimal flight route, which does not require huge 
computation times or a supercomputer to perform the analy- 
sis. The NOC analysis uses a feed-forward and linear 
feedback procedure for computing a near-optimum route at 
the wind velocity vector changes along the route. The 
computations are performed for dimensionless variables so 
that a given NOC solution can be applied to many other 
problems through appropriate scaling. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIGS. 1, 3 and 5 graphically illustrate optimal flight 
routes computed using an NOC analysis for two represen- 
tative wind patterns, as a function of flight distance remain- 
ing. 

FIG. 2 is a block diagram illustrating one embodiment of 
the control algorithm used in the invention. 

FIG. 4 graphically illustrates some feedback gain values 
10 computed for the wind pattern of FIG. 3, for different 
fractions of distance remaining. 

FIG. 6 graphically illustrates another embodiment, in 
which an optimal solution to an intermediate waypoint is 
extended to a subsequent waypoint. 

DESCRIPTION OF BEST MODES OF THE 
INVENTION 

FIG. 1 illustrates a windy environment in which an 
aircraft 11 operates. The aircraft 11 has a velocity vector 
V„=(V*cos 0, V sin 0) in a plane with Cartesian coordinates 
(x,y) and is subjected to a wind velocity vector field V w = 
(u(x,y),v(x,y)). The aircraft airspeed IVI=V, assumed con- 
stant here, is measured relative to the air mass, and 0 is a 
heading angle, measured in a plane and used as a control 
variable to minimize a time required to travel from a first 
waypoint (source) to a second waypoint (destination). 

Neighboring Optimal Control (‘NOC”) is a time varying 
linear feedback control approach that begins with a nominal 
30 solution (correct in the absence of other perturbations) and 
minimizes a performance index to second order in certain 
perturbed variables and parameters. In an NOC approach, 
feedback gains are variable and are tabulated as functions of 
a time-to-go parameter i (or of an along-track distance-to-go 
35 parameter). A perturbation in accumulated time from source 
to destination is estimated using a performance index J. The 
formalism is discussed by A. E. Bryson and Y. C. Ho, in 
Applied Optimal Control : Optimization, Estimation and 
Control Hemisphere Publ. Corp., New York, 1975, pp. 
40 148-158, 179-203. The invention extends this formalism by 
modeling winds as additional states of the system. 

FIG. 2 illustrates a feedback/feedforward configuration 
used to control a state vector perturbation 8x=<Sx, 8y, etc.) 
in the presence of variable cross-track wind shear, using an 
45 NOC approach. The control vector perturbation in a simple 
situation is 8u=(80), where 0 is a heading angle 0 for the 
aircraft A nominal solution (flight route, aircraft heading 
angle, etc.) is provided by a nominal solution module 21, 
with associated state variables x ru>m and wind velocity 
50 variables (u,v)^ m as a function of a route variable, such as 
time or time-to-go. The actual or assumed wind pattern 
(u(x,y),v(x,y)) and other relevant control variables are pro- 
vided by a dynamics module 29. Nominal solution state 
vector values and control vector values are computed and 
55 provided by the module 21 and are received at a first input 
terminal (+) of a difference module 23 (state values) and at 
a first input terminal (+) of a sum module 27 (control 
values). The dynamics module 29 generates and issues a 
state vector x that is received at a second input terminal (-) 
60 of the difference module 23. The first difference module 23 
forms and issues a difference state vector Ax=x nom -x that is 
received by an NOC gain matrix module 25. The gain matrix 
module 25 multiplies pre-computed NOC gains by compo- 
nent values of a state difference vector Ax and generates and 
65 issues a control difference vector Au, which is received at a 
second input terminal (+) of the sum module 27. The sum 
module 27 forms and issues a (corrected) control vector 
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u=u «>m + Au, which is received by the dynamics module 29 
and used to determine a new output signal, the state vector 
x. 

The system is initially expressed in spherical coordinates 
(0,d>) on a (nearly) spherical Earth. One then rotates the 
coordinate system to a new system (O mr <P mt ) to which a 
nominal path great circle route lies on the equator (0 TOf =O) 
in the new “normalized” coordinate system and the route 
destination corresponds to the azimuthal or longitude angle 
d> TO/ =0. For a route from point B (0 B & B ) to point A 
(0 A ,d> A ), one first computes an intermediate angle 

tan 2{-cos © A sin ©,+sin © A cos B g cos(<l> A -4>*), sin(<I>*- 


<J> A ) cos (1A) 

Intermediate dimensionless variables 

4i=sin(O A -O,.)c06 (IB) 

$ 2 ^cos 0^ sin 0/,-sin cos 0 F cos (1C) 

cos P+<t >2 sin P), (ID) 

> t>i c° s $-$2 s ^ n PX (IE) 


sin 0/.+COS Q A cos 0^ cos^-O,,)}, (IF) 

are then computed, where R is the distance (assumed to be 
constant) from the center of the Earth to die aircraft location 
(point P) in flight The latitude and longitude (angular 
variables) (B^.Op) to a general point become 

0^=-sin -1 (y' M X (1G) 

<jy^ tan 2 {S^-z’kY (1H) 

Dropping the subscripts on the coordinates (0^,d>^), 

equations of motion for the aircraft with heading angle 6 and 
wind velocity vector (u,v) become 


COS0 + tj ¥ (8) 

/ = sinfl + rjy 
5 0, 


where f is a dynamics vector depending upon the values of 
the state vector x f =(x\y’) and values of the control vector u 
10 (including the heading angle 0) and “0” represents one or 
more zero components of the dynamics vector f, such as 
wind vector component values at selected intermediate 
locations along the route. 

In determining an optimal solution in the presence of 
15 variable wind, the (transformed) path is decomposed into a 
grid of end-to-end segments, and the normalized wind 
velocity vector components, q M and are assumed to be 

piecewise linear within each segment, (i— 1)/(N^,— l)^x*^i/ 
20 (N*,-l) (N w =N tt , , N>, or N Kf ) for i=l, . . . , N^-l, between 

two consecutive grid points: 

y+n»(* x (9) 

n ,/*',y )=-n „/*') y+n (io) 

25 

Within this range of the coordinate x\ a piecewise linear 
wind velocity component q^x’.y 1 ) (w=u, us, v, vs) can be 
expressed as 

30 (ii) 

where i and q„, /+1 are the values of the variable q at grid 
points x’Ki-D/^-l), x’^i/fN^-l), and N u and N v are 
optionally different integers. 

35 Equations (5) and (6) become 

dx'/dt=cos 0-T] Mt (y)-y'+‘n M (j('), (12) 


dQ/dt={ V cos Of «}/(/? cos O}, 

(2) 


dy/0/*=sin MwCOZ+TUA 

(13) 

d<b/dt={ V sin 0fv}/R. 

(3) 

40 

The dynamic model is augmented with constraints 


Equations (2) and (3) can be simplified by using the 



(14 A) 

approximations cos d>«l and introducing the definitions 



(14B) 

x~R* p. 

(4A) 

45 

dr(^0^=O, 

(14C) 

y-/?e, 

(4B) 


to the values at each grid point. The first and 

second 




derivatives of the system equations are now piecewise linear 

L/R=0 RKf =mitvi\ longitude in rotated coordinates, 

(4C) 


functions of the coordinate y\ Other wind models might be 
used. Because of the linearization bv NOC, only the linear 

(4D) 

50 

variations remain. 


y'^y/L 

(4E) 


The system has a corresponding Hamiltonian 



(4F) 


H=XjJc os Of TL,(* y)}fXy(s in 0 +r|y(jc',y))fl, 

(15) 


(4G) 

55 

where A* and X y serve as Lagrange multipliers, to be 
determined, and q M =i^=0. The Euler-Lagrange equations of 

T[ v =v/V, 

(4H) 


motion for the system (H. Goldstein, Classical Mechanics , 
Addison-Wesley Publ., 1953, pp. 36-38, 215-234) to this 

Equations (2) and (3) become 



system become 


dx'/dt=cos 0f-T\,, 

(5) 

60 

dXJdt , =-bH/dx , ~-X^\\Jdx'~X^r[^x', 

(16) 

dy'/dt=sm 0 +r| w 

(6) 


dXy/dt’^—dH/dy '~~XJ^ l /dy'—X y dx\ Jdy 

(17) 

ox/di=j\x', a /*), 

0) 


0 =dHldQ=(—X x sin 0f cos 0). 

(18) 

The nominal path and nominal solution correspond to 0=0, 

65 

Equation (18) requires that 


y’=o. n u =nv=o> 



X y Jcos Q-Xjs'm 0=0. 

(19) 



US 6,600,991 B1 


5 

The Hamiltonian H is not an explicit function of time, and 
H=constant is an integral of the system. Here, the constant 
is set equal to 0. Equations (15) and (19) have the simulta- 
neous solutions 

Xv=-oos 6/{ l+n« cos 0+f|, sin 0}, (20) 


6 

specific values of the (dimensionless) coordinate x\ These 
perturbation values are then tabulated for look-up during 
application of the NOC approach. The final expression for 
computing a perturbed optimal heading angle 6 from the 
5 feedback gains and state perturbations is expressed as a 
single component vector 


Ay=— sin 0/{ 1+Tk cos 0+Tj, sin 0}, 

(21) 

One first computes a solution of a nominal trajectory, 
corresponding to the following choices: 

vii 

VII 

O 

(22 A) 


(22B) 

n^./*y) = *X M ' = H v > i* • • • » w w -i). 

(22C) 

0^=71. 

(22D) 


For the nominal solution of Eqs. (22AH22D), two of the 
Lagrange multipliers are found to be constants: 

(23A) 

(23B) 

A vector L(t') of other Lagrange multipliers at the final time 
t'=t'f satisfies additional boundary conditions: 

(24) 

where vp is a column vector representing final state 
constraints, v is a column vector representing undetermined 
multipliers, x' is a (dimensionless) vector including the 
location and wind parameter states, t'f denotes time corre- 
sponding to the final state, the superscript ‘T” indicates 
vector transposition, and 

^ = jc / -y/ (25) 

/-//. 

The partial derivatives of the state vector components are 
easily computed, and the Lagrange multipliers have the 
following values for the nominal solution: 



(26A) 

Vy=Ay= 0, 

(26B) 

A^(l SiSAU 

(26C) 

X vi =0(lSi£JV r ). 

(26D) 



SS=-K iu &x', 

(27 A) 

10 


(27B) 


where the feedback gain matrix coefficients, K^. 
are given generally by 

and K^, 



(28 A) 

15 

K ry =<m T lay-(n T la)Q-'R T , 

(28B) 


The quantity H^ 1 is generally a matrix (here, a scalar) 

20 

whose entries are the second derivatives of the Hamiltonian 
H in Eq. (15) with respect to components of the control 
vector u (here, the scalar heading angle 0), f tt is a vector 
whose entries are the partial derivatives of the dynamics 
vector f with respect to components of the control vector u 


(here, the heading angle 0). 


25 

The matrices Q A , R A and S~ are computed according to 


Q-=Q-(n T n)la, 

(29A) 


R=R-(m T n)la, 

(29B) 

30 

S'=S-(m T mya, 

(29C) 


and Q, R, S, m, n and a are matrix solutions of the following 


equations: 


35 


(30) 



(31) 



(32) 

40 

Cl(x\ u, H t'=ff)=(d®/df+L) f:/ =0. 

(33) 



(34A) 


dS/df — SA-A T S+SBSr-C, 

(34B) 

45 


(34C) 


dR/dt=-{A T -SB)R, 

(35A) 


R(Sj)=(dy/dx) 0 

(35B) 

50 

dQ/dt'-R r BR, 

(36A) 


A Backward Sweep method, discussed by A. E. Bryson 
and Y. C. Ho, op cit, is used here to obtain an NOC solution. 
Other methods that can be used to compute feedback gains 
and obtain nominal optimal solutions include Newton, 
Newton-Raphson, quasi-Newton, gradient algorithms and 
shooting algorithms. These methods can be used to obtain an 
optimal solution for nominal conditions (e.g„ with zero 
wind or a constant wind vector). Each state parameter is then 
perturbed slightly while holding each of the other param- 
eters constant. By computing the resulting changes in 0 and 
in tf corresponding to changes in the state parameters (e.g., 
positions and wind vectors), the neighboring optimal feed- 
back gains may be computed numerically, using a pertur- 
bation in 0 or tf divided by a perturbation in each of the state 
parameters. 

One result of application of the Backward Sweep method 
is numerical values of the perturbation feedback gains at 



e<tf)= o. 

(36B) 



(37 A) 

55 

m(t'f)^=(Bn/dx) r /, 

(37B) 


dn/dt’=R r Bm, 

(38 A) 


n(t'fyddy!dt), f 

(38B) 

60 

da/dt'=m T Bm, 

(39A) 



(39B) 


Development and use of these relations are discussed in 
Bryson and Ho, op dt, pp. 179-203. 

65 Feedback gain coefficients Ky(x’), K*, K^x’), 

(x') and K^x') are computed for a sequence of selected 
values of the nominal variable x'(0=x’^l), such as x'=0. 
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0.005, 0.01, 0.015, . . . , 1.000. and are provided as lookup 
table values. Any reasonable grid size will suffice here, 
provided the quantization errors are kept within tolerable 
limits. The feedback relation (27) is used here to compute a 
perturbation in the heading angle (control variable) 0: 

N (40) 

Afl = -Kj M - Ky’Uy' - £ K„j&{du./d?) l - £ tf^Av),. 

tel tel 

NOC routing with variable winds uses linear feedback of 
perturbations in the wind field (u,v) along a flight path to 
modify the aircraft heading angle to achieve a minimum 
time flight path. The algorithm uses the following procedure. 

(1) Convert aircraft latitude and longitude into rotated 
angular coordinates so that a nominal flight path lies on a 
great circle on the equator, with the destination waypoint 
corresponding to the rotated prime meridian. 

(2) Introduce dimensionless, normalized location coordi- 
nates (x',y’), with O^x’^l and y’=0 corresponding to a 
nominal solution flight path (optional). 

(3) Compute state perturbations for the present aircraft 
and wind states, vis-a-vis the corresponding nominal solu- 
tion values. 

(4) Interpolate feedback gain coefficients K from lookup 
tables. 

(5) Multiply the gain coefficients by the corresponding 
state perturbations and sum the products, as indicated in Eq. 
(40). 

(6) Add the heading angle perturbation to the nominal 
heading angle to obtain a resulting heading angle in the 
dimensionless coordinate system. 

(7) Convert to a resulting heading angle in the original 
(dimensioned) coordinate system. 

(8) Repeat the steps (3)-{7) for selected points along the 
route until the destination waypoint is readied. 

Step (2) is optional, and the analysis can be carried out in 
standard (non-normalized) coordinates. Equations (1A )- 
(1H) and (4A)-(4E) are applied to convert to normalized, 
dimensionless coordinates. With reference to step (3), the 
primary perturbations relative to a nominal solution path are 


A/=y'-0(i»m), (41 A) 

AT Ur=TV/*=l» • • • , AU (4 IB) 

, AT,)- (4 1C) 


The wind components typically are provided at selected grid 
points and are resolved along North and East directions and 
must be rotated and translated in order to express the wind 
components in a rotated coordinate system, such as the 
normalized coordinate system 
FIG. 1 illustrates a wind-optimal trajectory that might be 
experienced by an aircraft flying from the East Coast to the 
West Coast of the United States, or inversely, with flight 
parameters 


V=480 knots (552.4 miles per hour). 

(42A) 


(42B) 

v(*;y)=o, 

(42C) 

VA=« J =0.2(A=2400 nmi). 

(42D) 


In this example, only one measurement point is provided for 
u s . In this implementation of the invention, the nominal 
solution corresponds to u,=0, and the perturbation is a 
constant, non-zero wind shear parameter us, often measured 


8 

or retrieved just before a flight begins. For V=480 knots, 
U/=0.2 knots per nautical mile (nmi), and a straight line 
flight path distance of 2500 nmi, the reduction in flight time 
over the flight path is 11.8 minutes or 4 percent, relative to 
5 the time required for a constant heading flight (300.6, min 
versus 312.4 min). FIG. 1 graphically illustrates the constant 
heading flight path and an optimal solution flight path for 
this situation. 

In a second example applying the NOC approach, 
10 expressed in rotated, dimensionless coordinates and illus- 
trated in FIG. 3, along-track and cross-track wind velocities 
are both non-zero: 



£=2500 nmi. 

(43A) 

15 

V=480 knots, 

(43B) 


«(*, >) = (43.6 - 0.4y) kn (* = £) 
= 43.6 kn (* = £/2) 

(43C) 

20 

= (43.6+ 0.4>>) kn(* = 0). 



vt*iy)=-52.9 kn. 

(43D) 


25 where the wind velocity components, u and v, are interpo- 
lated linearly in each of the regions 0<x<L/2 and L/2<x<L. 
Table 1 and FIG. 3 compare the values and the routes for a 
nominal solution flight path and for an optimal solution 
flight path using the NOC approach. Here, flight along a 
30 nominal solution route (y'=0) requires 346 min, and flight 
along an NOC-optimized route requires 341.1 min, a dif- 
ference of 4.9 min or about 1.4 percent 

TABLE 1 

35 

Nominal and perturbed values of wind parameters 


«wi V W4 



Vw.,, 

(kn/nmi) 

V.. 3 

(kn/nmi) 

V..4 

(kn/nmi) 

(i=l, 2, 3) 
(kn) 

(i = 1, 2, 3) 
(kn) 

Nom- 

inal 

0 

0 

0 

39.9 

12.5 

Per- 

turbed 

0.4 

0 

-0.4 

43.6 

-52.9 


45 FIG. 4 graphically illustrates variation of a gain matrix 
coefficient K,, , with time-to-go parameter for three locations 
along a flight route, corresponding to source (i=l), mid-way 
(i=2) and destination (i=3). Each of the coefficients goes to 
0 as the time-to-go parameter decreases to 0. The feedback 
50 gain corresponding to a perturbation in wind velocity at or 
near the destination clearly affects the entire flight route, 
where total flight time is to be minimized. 

In a third example applying the NOC approach, expressed 
in standard coordinates and illustrated in FIG. 5, along-track 
55 and cross-track wind velocities are similar to the second 
example, but the flight distance is longer: 


£=6000 nmi, 

(44A) 

V=480 knots. 

(44B) 

u(x, y) - (43.6 - O.ly) kn (* = £) 

(44C) 

= 43.6 kn (x = L/2) 


= (43.6+ O.ly) kn (jc = 0), 


v(*y)=- 52.9 kn. 

(44D) 
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where the wind velocity components, u and v, are again 
interpolated linearly in each of the regions 0<x<L/2 and 
L/2<x<L. Flight along a nominal solution route (y-0) 
requires 830.5 min, and flight along an NOC-optimized 
route requires 803.3 min, a difference of 27.2 min or about 5 
3.4 percent 

In the second and third examples, because the NOC gains 
are computed in dimensionless form, a family of NOC 
solutions can be generated or approximated by scaling, 
using the scaling parameters L, V and V w , among others. As 10 
an example of such scaling, the solution in the second 
example can be reconfigured to provide an optimal solution 
for a flight over 6000 nmi (third example) between source 
and destination waypoints with V=480 knots and V ws =±0.1 
knots/nmi. The dimensionless feedback gain coefficients 15 
remain the same, but the perturbation values are scaled 
according to Eqs. (4A)-(4E). The dimensionless shear per- 
turbation at the origin is changed from (0.4)(2500)/480= 
2.0833 to (0.1)(6,000)/480=1.25. In this mapping of the 
second example parameters onto the third example 20 
parameters, the flight time difference (nominal solution 
flight path versus optimal solution flight path) changes from 
4.9 min (second example) to 27.2 min (third example). 

The prescription of zero wind velocity for a nominal 
optimal solution and neighboring optimal gain computation 25 
can be replaced by another wind profile, such as constant 
wind velocity, constant cross-track wind shear or another 
generalization. 

Where constant wind velocity is prescribed, the nominal 
optimal solution is easily computed. The heading angle Q nam 30 
is perturbed slightly to cause the aircraft to maintain a great 
circle path. If the along-track and cross-track wind velocities 
are uO and vO, respectively, the nominal heading angle in the 
rotated coordinate frame is given by 

35 

(45) 

where V is the (constant) aircraft speed, and the effective 
ground speed becomes 

V^VcosS^+uO. (46) 40 

One may use a more general wind profile when computing 
the nominal optimal solution, using any desired numerical 
optimization algorithm, such as a gradient method or 
Newton-Raphson method. For this more general profile, an 45 
analytical solution may not be possible, but the numerical 
solution is as straightforward as the solution for zero wind 
velocity. 

An advantage in using a wind profile other than the zero 
wind prescription is that the neighboring optimal solution 50 
may be mere accurate. A disadvantage in using a wind 
profile other than zero wind prescription is that the nominal 
optimal solution and corresponding neighboring optimal 
feedback gains must be recomputed for each specific way- 
point pair and wind condition. 55 

One may apply iterations to a neighboring optimal aircraft 
routing (NOAR) algorithm Here, one assumes zero wind 
velocity to compute the nominal optimal control and the 
corresponding neighboring optimal feedback gains. The 
neighboring optimal solution is then computed as indicated 60 
in the preceding development by incorporating the actual 
wind conditions measured for the particular flight The 
computed neighboring optimal trajectory is used as a (new) 
nominal optimal solution, and the NOAR algorithm is 
applied again, but with the perturbations in wind and wind 65 
shear computed with respect to the (new) nominal optimal 
solution path. The resulting heading angle perturbations are 
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now added to the nominal heading computed during the 
preceding iteration. This iterative process is repeated until 
the solution has converged to within a selected maximum 
inaccuracy. In many instances, a first iteration provides a 
(new) solution that is already sufficiently accurate, although 
additional iterations can be applied. 

Another extension to the NOAR algorithm restarts the 
NOAR computation after each iteration during the heading 
perturbation calculations. Beginning with the aircraft at the 
source waypoint (x=l, y’=0), from a time t'=0 (position 0 in 
FIG. 6) to a time t=At'(>0), one computes and applies the 
neighboring optimal heading perturbations as described in 
the preceding to bring the aircraft to a new location with 
coordinates x’=l-Ax’ and y’=Ay’ (position 1 in FIG. 6). A 
new coordinate rotation and normalization is performed to 
transform to a new normalized (dimensionless) coordinate 
frame for which x’ /MW =l and y' nw =0. The normalized time 
is also reset to t’^M). This process is repeated to move the 
aircraft to a new resulting location (position 2 in FIG. 6). The 
number of iterations is not limited to 3. This approach allows 
substantially exact feedback gains to be computed at each 
iterative step, rather than being obtained from a lookup 
table. This approach effectively reduces the scale of the wind 
grid as the aircraft approaches the destination waypoint so 
that the wind grid resolution improves as the aircraft 
approaches the destination. 

The aircraft speed V need not be constant throughout the 
flight. If one uses a variable speed profile, the aircraft speed 
to be applied at a given location (x,y) must be used in the 
coordinate normalization process. With that minor modifi- 
cation incorporated, the remainder of the NOAR algorithm 
is applied as discussed in the preceding. For example, if an 
aircraft speed Vj applies at a coordinate value x y <j=l, 2, 3), 
the wind components would be normalized as follows: 
r|y=vj/Vj(j— 1, 2, 3), where vj refers to a wind velocity 
component. 

The NOAR can be extended to include variations in wind 
velocity in the vertical coordinate. One difficulty with 
including the vertical coordinate is that commercial jet 
aircraft performance varies strongly with altitude. Although 
a primary goal of the NOAR algorithm is to determine a 
minimum-time path between the source and destination 
waypoints, the overriding goal is to determine a minimum- 
fuel-usage path. Fuel minimization is assumed to have been 
performed by choosing an appropriate flight altitude and 
airspeed V. By performing a minimum-time optimization for 
the selected fixed altitude, fuel consumption is further 
reduced. 

A straightforward extension would include aircraft alti- 
tude as a control variable and to introduce a three- 
dimensional wind model. The resulting vertical wind profile 
would be evaluated in terms of efficiency of aircraft and 
engine performance. Another level of complexity in this 
extension would include state models for the aircraft fuel 
flow rate as a function of altitude and airspeed, when 
computing the NOC gains. This is a non-trivial problem. 

Wind shear and bias wind parameters are determined from 
measurements taken at a sparse grid of points adjacent to one 
or more candidate flight routes. If the normalized wind 
velocity parameters rj u and r|„ (introduced in Eqs. (5) and 
( 6 )) are determined at two locations (xl',yT) and (x2\y2')4 
the corresponding wind shear parameters at an intermediate 
location (x3\y3’), with xl'^x , 3'^x2’ and yT^y3'^y2', can 
be estimated as 

i) tt> (x3 , ,y3 , )={n ii (jc',yr)-TL(Jf'. y2')}/(yl'-y2'), (47) 

y2')}l{y\'-y2'). (48) 
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Corresponding bias wind parameters can be estimated as 


(49) 

(so> 

Another approach for estimating wind shear and wind bias 
parameters fits two-dimensional polynomials, p tt (x’,y’) and 
p^x’y) for measurements, tl M (x\y') and t| v (x\y'), at several 
adjacent locations (x',y’>, estimates the wind bias values as 



(51) 10 

T\ v (x3\y3')=p v (x3\y3'), 

(52) 

and estimates the wind shear values as 


} (ja ^ 

(53) 15 


(54) 


The normalized coordinates values (x3\y3’) can be esti- 
mated using coordinates values for a preceding location on 
the flight path and a heading angle for that preceding 20 
location. Other interpolation methods can be used to esti- 
mate the wind bias and wind shear parameters from a sparse 
grid of measurements. 

The wind shear values, which are modeled as constants in 
Eqs. (9), ( 10) and 14A)-( 14C), can be generalized further as 25 
first order lags or as sinusoidally varying functions of the 
cross-track coordinate y\ However, use of constant wind 
shear values (zero or non-zero) is expected to be more useful 
in a general situation. 

The NOAR computations can be made on board the 30 
aircraft, before and/or during the flight, car can be made at 
another (“off-board”) computing facility not located on the 
aircraft. Alternatively, a combination of on-board and off- 
board computations can be made. This combination 
approach would be useful where updated or corrected wind 35 
information becomes available after the aircraft has 
departed. An off-board facility can receive updated wind 
information and can perform a selected preliminary portion 
of, or all, NOAR computations and can transmit Are raw data 
and/or the results of these computations to the target aircraft, 40 
for use in computing and adjusting the results of the modi- 
fied NOAR computations to provide an adjusted NOAR 
flight route; or the aircraft can obtain wind measurements 
locally and can provide some or all wind information 
updates itself as the flight progresses. 45 

The NOAR algorithm may be applied in at least two 
distinctly different modes. In the first of these modes, the 
calculation of a neighboring optimal trajectory in either 
on-board or off-board systems may be made just prior to 
flight so that the aircraft can be directed to follw the 50 
pre-computed trajectory. The neighboring optimal trajectory 
may be updated as often as desired as new wind information 
becomes available. 

In the second mode, instead of computing the entire 
neighboring optimal trajectory prior to flight in either 55 
on-board or off-board systems, one may use the neighboring 
optimal wind routing algorithm in a real-time feedback 
mode. In this mode, the computer would continuously 
calculate the current optimal heading to use to cause the 
aircraft to ultimately follow a neighboring optimal trajec- 60 
tory. 

Where relatively few intermediate waypoints are 
specified, aircraft personnel may be able to use a hand held 
specially programmed computer or one or more nomographs 
or dedicated slide rule mechanisms to determine heading 65 
angle adjustments and the like at selected locations along the 
flight route. 


What is claimed is: 

1. A method for determining an optimal route for an 
aircraft in moving between first and second waypoints in the 
presence of variable winds, the method comprising: 

providing a nominal solution for a route of an aircraft 
moving with a selected speed between the first and 
second waypoints in the presence of a first wind 
environment for which the nominal solution can be 
determined; 

providing values for a second wind environment that 
differs from the first wind environment; and 

using at least one computer to determine a neighboring 
optimal control solution for an aircraft route between 
the first and second waypoints in the presence of the 
second wind environment 

2. The method of claim 1, further comprising expressing 
said first and second wind environments and said location 
coordinates for said nominal route and for said optimal route 
in dimensionless variables. 

3. The method of claim 1, further comprising providing 
values for wind velocity in said fist wind environment that 
are substantially constant. 

4. The method of claim 1, further comprising providing 
values for wind velocity in said first wind environment that 
are substantially zero. 

5. The method of claim 1, further comprising providing 
values for cross-track wind shear in said first wind environ- 
ment that is substantially constant 

6. The method of claim 1. further comprising providing 
values for wind velocity in said second wind environment 
varies linearly with at least one location coordinate. 

7. The method of claim 1, further comprising selecting 
said aircraft speed to be constant 

8. The method of claim 1, further comprising selecting 
said aircraft speed to have at least first and second distinct 
values at first and second locations alone said route. 

9. The method of claim 1, further comprising, providing 
values for wind velocity in said second wind environment 
varies linearly with at least one location coordinate within 
each of at least two segments that extend between said first 
and second waypoints. 

10. The method of claim 1, wherein said step of providing 
said neighboring optimal solution comprises: 

converting aircraft latitude and longitude coordinates into 
rotated angular coordinates so that said nominal solu- 
tion lies along a great circle for the Earth; 

computing perturbations in at least one variable describ- 
ing a state of said aircraft for said optimal solution 
relative to said nominal solution, and in at least one 
variable describing said second wind environment rela- 
tive to said first wind environment; 

estimating at least one feedback gain coefficient for the at 
least one variable for the state of said aircraft; 

multiplying the estimated at least one gain coefficient by 
the perturbation in the at least one variable for the state 
of said aircraft to obtain at least one product value; 

summing the product values for the at least one variable 
for the state of said aircraft at at least one selected 
location along said route; 

using the sum of product values to obtain a heading angle 
perturbation of said aircraft for the at least one selected 
location: and 

combining the heading angle perturbation to said nominal 
solution heading angle to obtain an optimal heading 
angle at the at least one selected location. 
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11. The method of claim 10, further comprising express- 
ing said first and second wind environments and said loca- 
tion coordinates for said nominal route and for said optimal 
route in dimensionless variables. 

12. The method of claim 1, further comprising using at 
least one of (i) a computer located on board said aircraft and 
(ii) a computer that is not located on board said aircraft to 
determine said neighboring optimal solution. 

13. The method of claim 1, further comprising providing 
information on said second wind environment, after said 
aircraft has departed from said first waypoint, for at least one 
of (i) a computer located on board said aircraft and (ii) a 
computer that is not located on board said aircraft, and 
determining an updated neighboring optimal solution. 

14. The method of claim 1, further comprising providing 
at least one of said values for said second wind environment 
using measurements of a wind value made at at least two 
locations adjacent to a flight route given by said nominal 
solution. 

15. A method for determining an optimal route for an 
aircraft in moving between first and second waypoints in the 
presence of variable winds, the method comprising: 

providing a first nominal solution for a route of an aircraft 
moving with a selected speed between the first way- 
point and a selected third waypoint, intermediate 
between the first and second waypoints, in the presence 
of a first waypoint first wind environment for which a 
first waypoint nominal solution can be determined; 

providing values for a first waypoint second wind envi- 
ronment that differs from the first waypoint first wind 
environment; 

determining a neighboring optimal solution for an aircraft 
route between the first and third waypoints in the 
presence of the fist waypoint second wind environment; 

providing a second nominal solution for a route of an 
aircraft moving with a selected speed between the third 
waypoint and a selected fourth waypoint, intermediate 
between the third and second waypoints, in the pres- 
ence of a third waypoint first wind environment for 
which a second nominal solution can be determined; 

providing values for a third waypoint second wind envi- 
ronment that differs from the third waypoint first wind 
environment; and 

determining a neighboring optimal solution for an aircraft 
route between the third and fourth waypoints in the 
presence of the third waypoint second wind environ- 
ment 

16. The method of claim 15, further comprising selecting 
said fourth waypoint to substantially coincide with said 
second waypoint 

17. A system for determining an optimal route for an 
aircraft in moving between first and second waypoints in the 
presence of variable winds, the system comprising a com- 
puter that is programmed: 

to provide a nominal solution for a route of an aircraft 
moving with a selected speed between the first and 
second waypoints in the presence of a first wind 
environment for which the nominal solution can be 
determined; 

to provide values for a second wind environment that 
differs from the first wind environment; and 

to determine a neighboring optimal control solution for an 
aircraft route between the first and second waypoints in 
the presence of the second wind environment. 

18. The system of claim 17, wherein said computer is 
further programmed to express said first and second wind 
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environments and said location coordinates for said nominal 
route and said optimal route in dimensionless variables. 

19. The system of claim 17, wherein said computer is 
further programmed to provide values for wind velocity in 

5 said second wind environment that are substantially con- 
stant. 

20. The system of claim 17, wherein said computer is 
further programmed to provide values for wind velocity in 
said second wind environment that are substantially zero. 

10 21. The method of claim 17, wherein said computer is 

further programmed to provide values for cross-track wind 
shear that are substantially constant. 

22. The system of claim 17, wherein said computer is 
further programmed to provide values for wind velocity in 

15 said second wind environment that vary linearly with at least 
one location coordinate. 

23. The system of claim 17, wherein said computer is 
further programmed to select said aircraft speed to be 
constant. 

20 24. The system of claim 17, wherein said computer is 

further programmed to have at least first and second distinct 
values at first and second locations along said route. 

25. The system of claim 17, wherein said computer is 
further programmed to provide values for wind velocity in 

25 said second wind environment that vary linearly with at least 
one location coordinate within each of at least two segments 
that extend between said first and second waypoints. 

26. The system of claim 17, wherein said computer is 
further programmed: 

30 to convert aircraft latitude and longitude coordinates into 
rotated angular coordinates so that said nominal solu- 
tion lies along a great circle for the Earth; 

to compute perturbations in at least one variable describ- 
ing a state of said aircraft for said optimal solution 

35 relative to said nominal solution, and in at least one 
variable describing said second wind environment rela- 
tive to said first wind environment; 

to estimate at least one feedback gain coefficient for the at 

least one variable for the state of said aircraft; 

40 

to multiply the estimated at least one gain coefficient by 
the perturbation in the at least one variable for the state 
of mid aircraft to obtain at least one product value; 

to sum the product values for the at least one variable for 

45 the state of said aircraft at at least one selected location 
along said route; 

to use the sum of product values to obtain a heading angle 
perturbation of said aircraft for the at least one selected 
location; and 

50 to combine the heading angle perturbation with said 
nominal solution heading angle to obtain an optimal 
heading angle at the at least one selected location. 

27. The system of claim 26, wherein said computer is 
further programmed to express said first and second wind 

55 environments and said location coordinates for said nominal 
route and for said optimal route in dimensionless variables. 

28. The system of claim 17, wherein said programmed 
computer is located at at least one of (i) a location on board 
said aircraft and (ii) at a location not on board said aircraft, 

60 to determine said neighboring optimal solution. 

29. The system of claim 28, wherein said computer is 
further programmed to receive information on said second 
wind environment, after said aircraft has departed from said 
first waypoint, and to determine an updated neighboring 

65 optimal solution. 

30. The system of claim 17, wherein at least one of said 
values for said second wind environment is provided using 
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measurements of a wind value made at at least two locations 
adjacent to a flight route given by said nominal solution. 

31. A system for determining an optimal route for an 
aircraft in moving between first and second waypoints in the 
presence of variable winds, the method comprising a com- 5 
puter that is programmed: 

to provide a first nominal solution for a route of an aircraft 
moving with a selected speed between the first way- 
point and a selected third waypoint, intermediate 
between the first and second waypoints, in the presence 10 
of a first waypoint first wind environment for which a 
first waypoint nominal solution can be determined; 

to provide values for a first waypoint second wind envi- 
ronment that differs from the first waypoint first wind 
environment; 

to determine a neighboring optimal solution for an aircraft 
route between the first and third waypoints in the 
presence of the first waypoint second wind environ- 
ment; 
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to provide a second nominal solution for an aircraft route 
between the third waypoint and a selected fourth way- 
point intermediate between the third and second 
waypoints, in the presence of a third waypoint first 
wind environment for which a second nominal solution 
can be determined; 

to provide values for a third waypoint second wind 
environment that differs from the third waypoint first 
wind environment; and 

to determine a neighboring optimal solution for an aircraft 
route between the third and fourth waypoints in the 
presence of the third waypoint second wind environ- 
ment 

15 32. The system of claim 31, wherein said computer is 

further programmed to select said fourth waypoint to sub- 
stantially coincide with said second waypoint. 

* * * * * 
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